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Abstract 

We have investigated the ground state structures of neutral and charged XH4( X=C, Si, Ge, 
Sn and Pb ) molecules using the first-principles electronic structure methods. The structure of 
positively charged molecules for X = Si, Ge, Sn and Pb is characterized by a severe distortion 
from tetrahedral structure and an unusual H-H bond while the negatively charged molecules get 
distorted by pushing two hydrogen atoms away from each other. However, CH4+ and CH4~ are 
exceptions to this behavior. We provide an insight into the symmetry breaking mechanism and 
unusual H-H bonding using simple electrostatic arguments based on the unequal charge distribu- 
tion on H atoms. Those charged molecules having unequal charge distribution on H atoms get 
distorted due to different electrostatic forces between the atoms. We show that the directional- 
ity and occupation of the highest occupied molecular orbital(HOMO) play an important role in 
creating charge asymmetry in these molecules. 

PACS numbers: 73.22.-f, 36.40.Mr, 36.40.Wa,36.40.Qv 
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I. INTRODUCTION 



Although neutral XH4( X= C, Si, Ge, Sn and Pb ) molecules have tetrahedral symmetry, 
charging these molecules breaks this symmetry. The symmetry breaking in general is associ- 
ated with a large structural distortion and in many cases unusual H-H bonding arrangement. 
Some of these clusters have been investigated earlier and the symmetry breaking has been 
attributed to the Jahn- Teller effect. In the Jahn- Teller effect IjJ, , electrons occupying 
a degenerate level favor a low symmetric structure by lifting the degeneracy. For exam- 
ple, the highest occupied molecular orbital(HOMO) of tetrahedral SiH4 is triply degenerate. 
When an electron is removed from SiH4, its structure gets distorted drastically due to the 
Jahn-Teller effect fl fl E | . 

The structure of SiH4+ has been investigated by a large number of workersj^, 
and it took almost a decade to settle its ground state structure. Now it is generally agreed 
that the ground state structure of SiHA^has C, symmetry Similarly, the structure of CH4"'" 
has been studied extensively Q Q, Q, Q Q, llq . though there is still ambiguity 



regarding its ground state structure. The calculations done by the density functional theorv( 
DFT ) methods show that the ground state structure of 0134"*" has symmetry ^, , 
while the Hartree-Fockf HF ) and post HF calculations show the ground state structure as 

nn 

C2i) |9|, ll6| . The lowering of tetrahedral symmetry of CH4+ and SiH4"'' has been attributed to 
the Jahn-Teller effect. However, the structure of CH4''" is quite different from that of SiH4+. 
It is difficult to understand the difference between the CII4"'" and SiH4"'" distortions only from 
the Jahn-Teller effect. This is because the Jahn-Teller theorem only suggests the existence 
of distortion and says nothing regarding the nature and magnitude of the distortion. 

In this paper we focus on the ground state structures of charged XII4 molecules using 
the first-principles electronic structure methods and provide an insight into the symmetry 
breaking mechanism. We find that upon charging most of the molecules undergo large 
structural distortion which can be attributed to the Jahn-Teller effect. Further, the nature 
of distortions of positively charged molecules is different from that of negatively charged 
molecules. Interestingly, most of the positively charged molecules have an unusual H-H bond 
and appear like transition states of unstable XH2"'" and H2 complex. We try to understand 
the distortions, unusual H-H bonding and stability of these molecules. Vibrational analysis 
shows that these molecules are not transition states. The calculated fragmentation energy 
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indicates that the stabihty of XH4"'" molecules decreases from CH4"'" to PbH4"'". 

Using various exchange-correlation functionals such as local density approximation(LDA), 
local spin density approximation(LSDA), gradient corrected approximation(GGA) and spin- 
polarized GGA(GGA-SP) we have investigated the effect of different occupations in charged 
tetrahedral XH4 molecules. Our calculations show that the configuration corresponding to 
integer occupancy and asymmetrical charge density leads to the lowest energy. We provide 
simple arguments in terms of charge asymmetry and electrostatic interaction between atoms 
which explain the nature of distortions in charged XH4 molecules. For example, we are able 
to explain why the nature of distortion of SiH4"^ is different from CH4"'' and why the positively 
charged molecules distort differently from the negatively charged molecules. We study the 
structural distortions of SiH4+ in detail using conjugate gradient method and show that the 
evolution of the structure is governed by the electrostatic interaction between the atoms. 
Our investigations reveal that the distribution of charges on H atoms plays an important role 
in distorting the charged molecules. If the charges on H atoms are not equal, the structure 
gets distorted due to different electrostatic forces on the atoms. We also discuss how the 
charge asymmetry is created in these molecules by the directionality and occupation of the 
HOMO. 

The plan of the paper is as follows. In section II we give computational details of the 
present work. Section III contains results and discussion and finally in section IV we give 
our conclusions. 



II. COMPUTATIONAL DETAILS 



We have investigated the structure of neutral and charged XII4 using Vienna ab-initio 
simulation package( VASP and Gaussian98 package. We have cross checked our 
calculations using both packages and find that the results obtained by both methods are 
in good agreement. A number of calculations were done using the unrestricted Hartree- 
Fock( UHF ) and DFT methods using the Gaussian98 package. In the DFT calculations we 
have used several exchange-correlation functionals such as Slater's exchange with the local 
correlation functional of Perdew( SPL ), Becke's exchange which includes gradient density 
with gradient corrected correlation functional( BPW91) and Becke's 3 parameter hybrid 
functional( B3LYP ) 3,Q- 

In these calculations we have employed 6-311g** basis set for 
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molecules having C and Si atoms and LanL2DZ basis set with effective core potential for 



other molecules. The plane wave ultrasoft 



pseudopotential calculations 



the local density approximation ( LDA ) 20|, |2l| employing a simple cubic supercell of 20A. 



17l | were done with 



The number of plane waves used in the calculation is set by energy cut-off ( Ecut). We 
have used 287 eV Ecut for neutral and charged CH4 molecules and 150 eV Ecut for other 
molecules. The use of supercell technique for charged systems causes electrostatic interac- 



tion within the images of supercells 
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23l |. Calculation of charged systems are efficiently 
handled in VASP by applying a background charge to maintain the charge neutrality and by 

n 

adding dipole and quadrupole moment corrections |17:]. The structure and ionization poten- 
tia. of S,H,^ predicted by us.ug th. procedure are irr Cpe.agreement wth those of earlrer 

calculations based on the HP, MP2 and LSDA methodsl 

III. RESULTS AND DISCUSSION 
A. Ground state structures 

Pig. 1 shows the ground state structures of XH4 molecules under neutral and charged 
states using the Gaussian98 package with BPW91 functional. We first discuss the ground 
state structures of the positively charged molecules. We see from the figure that the tetra- 
hedral symmetry is broken when XH4 molecules are ionized. The structure of CH4+ has 
symmetry, while other positively charged molecules have symmetry. We have investigated 
various possible symmetries of XH4"'" molecules and our results are summarized in Table. I. 
We see from the table that and structures of CH4"'" are very close in energy and differ 
only by 0.11 eV. The other positively charged molecules also have as the next higher 
energy structure. We could not find local minima of C-^v and symmetries for 0134"*" and 
C2V symmetry for SiH4+ at the BPW91 level. These missing structures of CH4+ and SiH4"'" 
molecules in our calculations have 
potentials and different methods [ 
state structures. 

It is interesting to note that in the symmetry, the positively charged molecules have 
a structure which is severely distorted from the tetrahedral structure. Purthermore, there is 
an unusual bonding between two hydrogen atoms. This bonding is confirmed by the charge 



been studied earlier using different exchange-correlation 
'P, I9, 15, 1^ which show that these are not the ground 
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density calculation as shown in Fig. 2 which shows large electron density between the two 
hydrogen atoms of SiH4+. To ascertain that this is not an artifact of the BPW91 calculation, 
we have re-optimized the ground state structures using the HF and DFT methods with 
different exchange correlation functionals such as SPL and B3LYP. The H-H bond lengths 
obtained from various methods are essentially the same. This kind of over-coordination of 
M.o.en ^ a>so been observed ,„ Cn^U. S.nM a.d ..a., .ydrogenated s.Ucon 
clustersj^]. In Fig. 3, we have shown the total energy curve of 81134+, H2 and H2"'" as a 
function of H-H distance. The total energy of SiH4"'" is calculated by keeping Si-H bond 
length fixed but varying H-H bond distance. From the figure we see that the curvature of 
SiH4"'" curve is smaller than H2 and larger than H2+ curves. This shows H-H bond in SiH4"'" 
is stronger than H2+ but weaker than H2. Vibrational analysis also indicates that the H-H 
bond is sufficiently strong in these molecules. We see from Table. II that the force constant 
for H-H stretch mode of H-H pair is smaller than H2 and larger than H2+ which implies that 
the strength of H-H bond is weaker than H2 and stronger than H2'''. 

For the ground state structures of the negatively charged molecules, we note from Fig. 
1 that the distortions of negatively charged molecules are different from those of positively 
charged molecules. XH4^ molecules have symmetry except CH4~ which has tetrahedral 
symmetry. The two hydrogen atoms of XH4~ molecules are pushed away from each other 
which is in contrast to the H-H bond of XH4+ molecules with Cs symmetry. 



B. Stability 

One may argue that the structure of XH4+ molecules with symmetry appears more 
like XH2+ and H2 complex and hence unstable. Frey and Davidson^] studied SiH4+ using 
CI singles root calculation and showed that structure is more stable than against 
fragmentation into SiH2"'" and H2. We have calculated the fragmentation energy which is 
the energy required to fragment a XH4"'" molecule into its binary products and has been 
calculated by taking the total energy difference between the XH4+ and the possible binary 
products. Of all the possible channels, the fragmentation of XH4~^ into XH2''' and H2 requires 
least energy. We give the fragmentation energies of XH4~*' molecules for this process in Table. 
III. We see from the table that the energy required for the fragmentation of XH4"^ into XH2"'" 
and H2 decreases for X=C to Pb. This means that the stability of XH4"'" molecules decreases 
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from C to Pb molecules. 

To examine the stability further we have done vibrational analysis of these molecules. 
In Table. IV we present the vibrational frequencies of charged and neutral XH4 molecules. 
Our calculated vibrational frequencies 880, 950, 2175 and 2192cm^^ of SiH4 are in good 
agreement with those of experimentally obtained valuesj3l 914, 953, 2189 and 2267cm~^. 
We see from the table that PbH4"'" has a negative frequency and other molecules have only 
positive frequencies. Further, we see that one frequency is very low for symmetric XH4"'" 
molecules and decreases from Sill4"'^ to GeH4"'" and become negative for PbH4"'". 

The positive frequencies of a molecule indicate that the structure corresponds to an 
energy minimum. If one frequency is negative, the structure may be a transition state of 
the molecule llBi]. Since the frequency calculation is highly dependent on the method, the 
small negative frequency of PbH4+ may be an artifact of the method. To clarify this we 
have done vibrational analysis on re-optimized structures of symmetric XH4"'" molecules 
using various methods. Calculation with the SPL exchange-correlation functional shows 
that one frequency of PbH4+ which was negative in the calculation with BPW91 functional 
is positive but near zero. Similarly, calculation with the UHF method also shows that all 
the frequencies of PbH4^ are positive. Irrespective of the sign of calculated frequencies, all 
these methods show that one vibrational mode of symmetric XH4"'" molecules has a low 
frequency. Also we find that the force constant corresponding to this mode is nearly zero, 
which implies that the H atoms in these mode are moving on nearly fiat potential energy 
surfaces. 

We see from Fig. 1 that the mirror symmetry plane of Cg structure is formed by the H-H 
pair and X atom. The vibrational mode with low frequency corresponds to the vibration of 
H-H pair about the mirror symmetry plane of structure. When the H-H pair vibrates 
about the plane, one atom of H-H pair moves in and the other moves out of the mirror 
plane without change in the H-H distance. Since the force constant is nearly zero for this 
vibration, the motion of H-H pair is more like an internal rotation about an axis which is 
perpendicular to the H-H bond and passing through the X atom. This internal rotation of 
H-H pair in XH4+ is similar to the internal rotation of 0-H in CH3OH molecule^, 30 



Thus the low frequency mode in the present case does not indicate a transition state. 
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C. Effect of different occupations 



Before proceeding further, we would like to resolve an interesting question regarding 
the occupation of degenerate states in charged tetrahedral XH4 clusters. Let us ionize a 
charged XH4 cluster keeping its tetrahedral structure intact. In this structure the HOMO 
is triply degenerate and is occupied by only 5 electrons. Hence, the HOMO of XH4"*" is not 
completely filled. In such case of incompletely filled degenerate systems, there are many 
possible ways to fill the levels. The question is what is the right way of filling the levels and 
which distribution leads to the lowest energy. The usual practice in such situation is to fill 
all levels with fractional number of electrons, thereby keeping the symmetry intact. 

To resolve the issue of occupation we have performed paramagnetic and spin-polarized 
calculations on tetrahedral SiH4+ and CH4"'" using VASP package with local and non-local 
functionals for various possible electronic distributions. In paramagnetic calculations we 
have chosen two possible electronic distributions. The configuration with integer ( 2,2,1 ) 
occupation corresponds to asymmetric charge density and the other with fractional( |, |, 
I ) occupation corresponds to symmetric charge density. Note that the spin multiplicity of 
SiH4"*' and CH4"'' is unity. The multiplicity is preserved in charge asymmetric configuration 
whereas it is zero for charge symmetric configuration. In spin-polarized calculations we 
have chosen three possible electronic configurations. The configuration with ||, ||, || 
( ||, ||, || ) occupation corresponds to charge symmetric density and multiplicity zero. 
Another configuration with 1|, 1|, 1|( ||, ||, ||) occupation also corresponds to symmetric 
charge density but multiphcity one. The other configuration with 1|, If, 1|( 1|, 1|, 0|) 
occupation corresponds to asymmetric charge density and multiplicity one. We note that in 
all these configurations, the Kohn-Sham wave function is a single slater determinant. Hence 
the energy of the multiplet can be calculated directly using von Barth scheme js^. The 
calculated total energy of tetrahedral SiH4+ and CH4+ molecules are summarized in Table. 
V. 

From the table we see that paramagnetic calculations with local( LDA ) and non-local( 
GGA ) functionals do not give much difference between total energies of the molecules for 
symmetric and asymmetric charge density. On the other hand, the spin polarized calcula- 
tions with local functional( LSD ) favors charge symmetry for SiH4+ and charge asymmetry 
for CH4~*'. However, for both the molecules LSD calculations show that the symmetric charge 
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density which has multiphcity of zero is highly unfavored. We note that in spin polarized 
calculations the exchange interaction is taken care of more accurately than in paramagnetic 
calculations. Hence the unphysical zero multiplicity is strictly not favored by spin polarized 
calculations. But LSD does not give consistent results for multiplicity preserved symmet- 
ric and asymmetric charge density of SiH4"'" and CH4+ molecules. This is a well known 



33l |. This deficiency of LSD 



shortcoming of the LSD in the context of atomic calculations 

is removed by the gradient corrected functionals 3] • We have therefore done spin polarized 
GGA( GGA-SP) calculations on SiH4"'" and CH4"'" molecules. Our calculations with GGA- 
SP functional give the consistent result that the asymmetric charge density is the ground 
state configuration for both SiH4"'" and CH4"'" molecules in tetrahedral structure. 

Apart from the total energy calculation, one more important issue is the calculation of 
forces in such systems. We find that enforcing charge symmetry with fractional |, |, | occu- 
pation in the degenerate level of XH4"'' in LDA calculations does not distort the tetrahedral 
structure. This is in contrast with the structural relaxation of tetrahedral XH4"^ with 2,2,1 
occupation of degenerate level where the molecules results in distorted structure. Although 
LDA does not distinguish clearly the charge symmetric and asymmetric configurations, the 
calculated forces are different for the two configurations. This indicates that by enforc- 
ing symmetric charge distribution, the forces on atoms of XH4"'" molecules are incorrectly 



estimated. It was also noted by Janak and Williams j33| that enforcing spherical charge 
distribution on atoms can give rise to errors in Hellman-Feynman forces. Therefore, our 
results clearly indicate that the electron charge density of tetrahedral XH4+ is asymmetric. 



D. Structural evolution 



To gain insight into the mechanism of structural distortion in SiH4+ we look at the evolu- 
tion of SiH4+ structure by using the VASP package with LDA functional and the conjugate 
gradient(CG) methodjs^. We first ionize SiH4 by holding atoms to obtain tetrahedral SiH4"'" 
and then relax the structure. We find that with 64 steps the system rolls down from the 
initial tetrahedral to a well converged ground state structure. Each CG step corresponds 
to a particular geometrical configuration. The bond length between silicon atom and the 
distance between the hydrogen atoms as a function of CG steps are shown in Fig. 4(a). Ini- 
tially one H atom moves away from the silicon atom and attains maximum Si-H bond length 
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of I.92A at 8th step. Then this repelled H atom moves closer to the nearest hydrogen from 
13th step which can be seen from the H-H distance shown in Fig. 4(b). The two hydrogen 
atoms come very close to each other to a distance of 0.84A and form an H-H pair. At this 
stage the structure of SiH4+ is C2V After 23rd step, this H-H pair gets rotated towards the 
other two hydrogen atoms without changing the H-H pair distance which finally results in 
Cs symmetric SiH4+. 

The summation of occupied Kohn-Sham orbital energies as a function of CG step is shown 
in the Fig. 4(c). From the figure we can see that the sum increases initially up to the 8th 
step and then decreases enormously between 13th and 23rd step. Then it fluctuates between 
24 and 40th step and eventually converges to a constant value. The figure shows that the 
structural changes occurring from 13 and 23 are due to lowering in the orbital energies. Note 
that there is a minimum at 23rd step. As mentioned earlier at this CG step the structure 
of SiH4"'" is with the H-H bond. In a sense, the sum of orbital energies shown in the 

n 

figure is similar to the Walsh diagram|36] where one takes into account only the sum over 
orbital energies in predicting the structure. In Walsh diagram, each orbital energy is plotted 
as function of bond angle of a molecule. Then the minimum of the sum helps in predicting 
the structure. Thus in the present case the Walsh diagram would predict the symmetry of 
SiH4"'' as symmetry and not C^, which is the correct symmetry. 



E. Role of electrostatic interaction 



We now show that the distortion and formation of H-H bond can be easily understood 
from arguments based on electrostatic interaction between atoms of the molecule. We first 
try to understand the structural evolution of SiH4+ as discussed in the previous subsection. 
Let us first consider neutral SiH4 molecule. Since H is more electronegative!^ than Si, the H 



atoms in SiH4 has a small negative charge w^ 
also shown by population analysis 



e charge 

30,0,1 



lile Si will have a small positive charge. This is 
I according to which each hydrogen atom holds 
a charge of -0.12e and silicon +0.52e. When the molecule is ionized to SiH4"'" holding its 
atoms fixed, population analysis shows that the hydrogen atoms have -0.12e, -0.02e, +0.06e 
and +0.30e charge and silicon +0.77e charge. Interestingly, hydrogen atoms of tetrahedral 
SiH4"'' do not have equal charge although the structure is symmetric. One hydrogen atom 
has more positive charge than the others and gets repelled by silicon atom which is also 



9 



positively charged. This hydrogen moves away from the sihcon to reduce the electrostatic 
repulsion and is attracted by the nearby hydrogen atom which has more electrons and forms 
an H-H bond. This H-H pair is positively charged and hence attracted by the nearby two 
hydrogen atoms which results in rotation of the pair towards the other two hydrogen atoms. 

Since the electronegativities of Ge, Sn and Pb atoms are lower than that of hydrogen 37| , 
GeH4+, SnH4+ and PbH4"'" have distortions and bonds similar to that of SiH4+. However, 
as noted earlier, the structural distortion of CH4+ is different from that of SiH4"'"; two 
hydrogen atoms of CH4+ move away from each other which is in contrast with the formation 
of H-H bond in SiH4^. We now show that the same arguments based on electronegativity 
and electrostatic interaction also explain the structural distortion of CH4+. Since carbon 
has larger electronegativity than H atom, carbon in CH4 will have negative charge while 
hydrogen will have positive charge. This is confirmed by population analysis which shows 
that a charge of -0.82e on carbon and +0.20e on each hydrogen atom. Ionizing CH4 in the 
tetrahedral symmetric structure results in +0.20e, +0.23, +0.42e and +0.55e charge on H 
atoms and -0.40e charge on carbon. We see that the charge distribution is asymmetric and 
two hydrogen atoms are more positively charged than others. The two H atoms repel each 
other and open H-C-H bond angle as seen in Fig. 1. 

Now we discuss the distortion of negatively charged molecules and explain it on the basis 
of the above arguments. We see from Fig. 1 that CH4~ has the tetrahedral symmetry 
while other negatively charged molecules have C2V symmetry. Putting an electron on CH4 
results in -l.OOe on carbon and nearly zero electronic charge on H atoms. Because of this 
symmetrical charge distribution on hydrogen atoms, the forces on H atoms are radial and 
the structure of CH4^ remain tetrahedral. In contrast, in other molecules the charge on H 
atoms are not equal. For example, in SiH4~ population analysis shows a charge distribution 
of -0.14e, -0.15e, -0.16e and -0.37e on H atoms while silicon holds a charge of -0.18e. We can 
see that the H atom which is more negative gets repelled by other atoms. As a result, one 
of the H-Si-H bond angle opens up with the corresponding expansions in Si-H bond lengths. 
The trend in the electron distribution of GeH4~, SnH4~ and PbH4~ is similar to that of 
SiH4~. Therefore, the distortion of these molecules are similar to that of SiH4~ and can be 
explained using the same reasoning. 
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F. Charge asymmetry 



It is interesting to note that in most of the charged molecules the symmetry is broken. 
We find that the symmetry breaking is associated with the creation of asymmetry in the 
charge distribution on H atoms upon charging. Removing or putting an electron on these 
molecules without changing the tetrahedral geometry may result in charge imbalance on 
H atoms or equal charges on H atoms as in the case of CH4^. This charge imbalance or 
balance on H atoms depends on the electronic structure of the molecule. In particular, on 
the nature and occupation of the HOMO. 



Using molecular orbital theory 



39| we discuss the nature of HOMO in these molecules and 



show how it contributes to charge distribution among H atoms. In Fig. 5, we show the MO 
construction of XH4 from linear combinations of atomic orbitals. As shown in the figure, 
there are eight MOs, the lower four are bonding states and the upper four are antibonding 
states. The bonding and antibonding states have either Ai symmetry or T2 symmetry. Ai 
level is non-degenerate while T2 level is triply degenerate. While the bonding energy level 
of Ai symmetry is always lower in energy than T2, the order of antibonding levels depends 
on the molecule. For CH4 antibonding Ai level is lower in energy than T2 while for SiH4, 
GeH4, SnH4 and PbH4 antibonding T2 level is lower than Ai. From the figure we can see 
that the MO with Ai symmetry is a linear combination of valance atomic s orbitals of X 
and H atoms. The MO with T2 symmetry is a linear combination of valance p orbitals of X 
and s orbitals of H atoms. Because of mixing between p and s character, T2 orbitals have 
strong directionality along H atoms. If T2 orbital is fully occupied, it has equal weight along 
H atoms and results in a uniform charge distribution on H atoms as in the case of neutral 
XH4 molecules. If T2 orbital is incompletely filled with unequal electrons, it has unequal 
weight along H atoms and results in a non-uniform charge distribution on H atoms. 

Removing an electron from tetrahedral XH4 results with 5 electrons and the T2 orbital 
is incompletely filled. Hence the charges on H atoms of tetrahedral XH4"'" molecules are not 
equal. This can be seen from Fig. 6(a) and (b) which show the charge densities of SiH4 and 
SiH4+ in tetrahedral structure. We see from the figure that while the charge densities on H 
atoms of SiH4 are same, the charge densities on H atoms of SiH4+ are different. In the case 
of negatively charged molecules the HOMO can be of T2 or Ai symmetry. The HOMO of 
CH4~ has Ai symmetry. Since Ai is a linear combination of valance s atomic orbitals with 
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equal weights from hydrogen atoms, the charge distribution is equal on H atoms of CH4~ 
as shown in Fig. 6(c). The higher electron density around the carbon atom as seen in the 
figure can be attributed to the higher electronegativity of carbon compared to hydrogen. 
The HOMO of other negatively charged molecules has T2 symmetry and is singly occupied. 
This leads to asymmetric charge distribution on H atoms. Note that since Ai level in CH4~ 
is non-degenerate, CH4~ does not show the Jahn- Teller effect, and also no charge asymmetry. 
On the other hand, if the HOMO has T2 symmetry and is not completely occupied as in the 
case of SiH4+, the molecule will show the Jahn-Teller effect and the charge asymmetry. This 
imphes that the creation of the charge asymmetry and the Jahn-Teller effect are intimately 
connected. 



IV. CONCLUSIONS 



We have investigated the ground state structures of neutral and charged XH4 molecules 
using first principle electronic structure methods. We find that charging XH4 molecules can 
break the symmetry and change their structures drastically. Negatively charged molecules 
are distorted by pushing two hydrogen atoms away from each other while the positively 
charged molecules get distorted with an unusual H-H bond formation except for X=C. The 
calculations done by various methods and the charge density analysis confirm the existence 
of H-H bond. Furthermore, vibrational analysis shows that the strength of H-H bond in the 
positively charged molecule is stronger than H2'*' and weaker than H2. We find that one 
vibrational mode of symmetric XH4+ molecules has low frequency and corresponds to 
internal rotation of the H-H pair. Our vibrational analysis indicates that the molecules are 
not transition states. The fragmentation behavior of XH4"'" into XH2 and H2 shows that the 
stability decreases from CH4+ to PbH4+. 

From the total energy calculations using LDA, LSD, GGA and GGA-SP functionals we 
have shown that the ground state electronic structure of a charged tetrahedral molecule 
has asymmetric charge density corresponding to integer occupation. We have provided an 
insight into the symmetry breaking and shown that the distortion and unusual H-H bond 
formation can be understood easily using simple arguments based on charge asymmetry and 
electrostatic interaction between the atoms. Unequal charges on hydrogen atoms give rise 
to different electrostatic forces on different atoms and lead to distortion. On the other hand 



12 



if the charges are equal on hydrogen atoms, the structural symmetry of the molecule will be 
preserved. We have shown how the charge asymmetry arises due to character of the HOMO 
and eventually leads to the symmetry breaking. 
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FIGURE CAPTIONS 



Fig. 1. Ground state structures of neutral and charged XH4 molecules. 

Fig. 2 Electron density contour of SiH4+ in the plane of Si and the two bonded H atoms. 

Fig. 3 The total energy of SiH4+, H2 and H2'^ molecules as a function of H-H distance. 

The total energy of the molecule is given with respect its ground state energy. 

Fig. 4. (a) The Si-H bond lengths and (b) H-H distances in SiH4"'" as a function of CG 

steps. The H atoms have the same label as those of SiH4"'" in Fig. 1. (c) The sum of energy 

levels of SiH4+ as a function of CG steps. 

Fig. 5. A schematic energy-level diagram showing the construction of molecular orbitals 
of XH4 molecule from atomic orbitals. The levels on the extreme left side are the valance 
states of X atoms and the level on the extreme right is the valance state of H atom. The 
levels which are in the center are the molecular states of XH4. 

Fig. 6. The isosurface contour of valance electrons of (a) SiH4 (b) SiH4"'" and (c) CH4~ 
obtained by using the Gaussian package. The value of the contour is chosen such that one 
able to distinguish the difference between the charge density on the H atoms. The densities 
shown in figure have the contour value of 0.20 for SiH4 and SiH4+ molecules and 0.3 for 
CH4-. 
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Table. I. Relative energies ( in eV) of XH4+ molecules for various symmetries with 
respect to the ground state structures. The results are obtained using the Gaussian98 with 
the BPW91 exchange-correlation functional. 







C>2d 


^3v 


c. 


CH4+ 


0.11 









SiH4+ 




0.73 


0.73 





GeH4+ 


0.83 


0.94 


0.87 





SnH4+ 


1.36 


1.46 


1.35 





PbH4+ 


1.80 


1.91 


1.72 






Table. II Properties of H-H bond calculated using the Gaussian98 package with the 
BPW91 exchange-correlation functional. 





SiH4+ 


GeH4+ 


SnH4+ 


PbH4+ 


H2 


H2+ 


H-H bond length( 1 ) 


0.80 


0.78 


0.77 


0.76 


0.75 


1.12 


H-H frequency ( cm""^ ) 


3598 


3911 


4061 


4134 


4349 


1976 


H-H force constant( eV/1) 


7.69 


9.08 


9.79 


10.15 


11.23 


2.32 



Table. III. The energy required to fragment XH4+ into XH2'^ and H2. The results are 
obtained using the Gaussian98 package with the BPW91 exchange-correlation functional. 



Reactant 


product 1 


product 2 


Fragmentation 

energy ( eV ) 


CH4+ 


CH2+ 


H2 


3.03 


SiH4+ 


SiH2+ 


H2 


0.69 


GeH4+ 


GeH2+ 


H2 


0.42 


SnH4+ 


SnH2+ 


H2 


0.28 


PbH4+ 


PbH2+ 


H2 


0.21 
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Table. IV The calculated vibrational frequencies u{ in cm~^ ) and force constant k( 
in eV/A ) of neutral and charged XH4 molecules using the Gaussian98 package with the 
BPW91 exchange-correlation potential. 





Neutral molecules 




CH4 


SiH4 


GeH4 


SnH4 


PbH4 










K 




1, 

K 




K 






1 


1298 


1.17 


879 


0.50 


788 


0.38 


676 


0.28 


644 


0.25 


2 


1298 


1.17 


880 


0.50 


788 


0.38 


676 


0.28 


644 


0.25 


3 


1298 


1.17 


881 


0.50 


788 


0.38 


676 


0.28 


644 


0.25 


4 


1521 


1.37 


949 


0.53 


901 


0.48 


735 


0.32 


703 


0.29 


5 


1521 


1.37 


950 


0.54 


901 


0.48 


735 


0.32 


703 


0.29 


6 


2975 


5.2(3 


2175 


2.81 


2058 


2.51 


185G 


2.05 


1769 


1.86 


7 


3090 


6.19 


2191 


2.97 


2071 


2.59 


1856 


2.07 


1786 


1.91 


8 


3090 


6.19 


2191 


2.97 


2072 


2.59 


1857 


2.07 


1786 


1.91 


9 


3090 


6.19 


2191 


2.97 


2072 


2.59 


1857 


2.07 


1787 


1.91 




Positively charged molecules 




CH4+ 


SiH4+ 


GeH4+ 


SnH4+ 


PbH4+ 




V 


k 


V 


k 


V 


k 


V 


k 


V 


k 


1 


470 


0.15 


256 


0.04 


136 


0.01 


46 


0.00 


-104 


0.01 


2 


471 


0.15 


616 


0.25 


510 


0.16 


380 


0.09 


316 


0.06 


3 


1022 


0.78 


680 


0.28 


532 


0.17 


394 


0.09 


339 


0.07 


4 


1253 


0.93 


766 


0.35 


542 


0.18 


430 


0.11 


396 


0.09 


5 


1420 


1.20 


850 


0.45 


658 


0.26 


495 


0.15 


425 


0.11 


6 


2G12 


4.26 


1077 


0.G9 


781 


0.37 


GGG 


0.27 


G31 


0.24 


7 


2758 


4.52 


2164 


2.83 


1957 


2.29 


1729 


1.78 


1506 


1.35 


8 


2851 


5.32 


2242 


3.15 


2046 


2.54 


1807 


1.96 


1623 


1.57 


9 


2852 


5.32 


3598 


7.69 


3911 


9.08 


4061 


9.79 


4134 


10.15 



18 



TABLE IV: Continued 





Negatively charged molecules 




CH4- 


SiH4- 


GeH4- 


SnH4 


PbH4- 




;/ 


k 


;v 


k 




k 


;v 


k 


// 


k 


1 


1133 


0.88 


690 


0.30 


614 


0.23 


491 


0.15 


347 


0.07 


2 


1133 


0.88 


759 


0.36 


729 


0.33 


619 


0.23 


469 


0.13 


3 


1133 


0.88 


790 


0.38 


768 


0.35 


645 


0.25 


673 


0.27 


4 


1372 


1.12 


955 


0.58 


884 


0.48 


784 


0.37 


675 


0.27 


5 


1372 


1.12 


1000 


0.59 


934 


0.52 


794 


0.37 


692 


0.29 


6 


2537 


3.82 


1323 


1.09 


1099 


0.73 


1067 


0.68 


766 


0.35 


7 


2631 


4.55 


1429 


1.22 


1323 


1.04 


1245 


0.92 


1124 


0.75 


8 


2631 


4.55 


1958 


2.32 


1786 


1.91 


1607 


1.54 


1472 


1.29 


9 


2631 


4.55 


1962 


2.37 


1804 


1.96 


1615 


1.56 


1489 


1.32 
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Table. V Relative energies( in eV ) of SiH4+ and CH4+ in tetrahedral structure for various 
electronic configuration with respect to charge asymmetric configuration, using different 
exchange-correlation functionals. 





configuration 


SiH4+ 


CH4+ 


LDA 


2, 2, 1 


0.000 


0.000 




5 5 5 
3 ' 3 ' 3 


0.001 


0.001 


GGA 


2, 2, 1 


0.000 


0.000 




5 5 5 
3 ' 3 ' 3 


0.030 


0.000 


LSD 


IT, IT, 1T( U, U, Oi) 


0.000 


0.000 




5-1- 5f 5-1-/51 5| 5|\ 
3 1 ' 3 1 ' 3 1 V 3-^' 3^' 3-^/ 


0.139 


0.296 




IT, IT, 1T( ii, fi, |i) 


-0.028 


0.036 


GGA-SP 


IT, IT, 1T( U, U, Oi) 


0.000 


0.000 




5-1- 5-|- 5-|-/5| 5|5|\ 
3 1 ' 3 1 ' 3 1 V 3^'3^) 


0.215 


0.333 




IT, IT, IT ( ii, ii, ii) 


0.070 


0.105 
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( PbUi ) ( PbHi+ ) ( PbH^- ) 




H-H distance (A) 



<^ 


1.9 


1 1 1 1 1 1 1 1 1 1 1 

A 


1 1 1 1 
Si-H4 

/ 


-(— > 


















1.7 




^ Si-H3 - 


o 






Si-Hl,Si-H2" 

r 


1 

1 


1.5 


-, \/\ /•■■■■■.. /■■/"'■■„•■ •■- 






i'...--' 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 










Z.'f 








1.8 


" V (b) 


^Hl-H2 


• 1— 1 






1 


1.2 


1 1 1 1 1 V 1 1^ i U- 


/H3-H4 ~_ 

— \ — \ — \ — \ — 




-104 


A Energy sum of occupied levels 


> 








^1 


-108 
















-112 









8 24 40 

CG steps 




X 



H 



